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' Abstract. Wc study the collision of two fast solitons for the nonlinear Schrodingcr 

O,. equation in the presence of a spatially adiabatic external potential. For a high 

<^ I initial relative speed ||ij|| of the solitons, we show that, up to times of order 

log ||w|| after the collision, the solitons preserve their shape (in L^-norm), and 

the dynamics of the centers of mass of the solitons is approximately determined 

by the external potential, plus error terms due to radiation damping and the 

(— I ' extended nature of the solitons. We remark on how to obtain longer time scales 

^ |. under stronger assumptions on the initial condition and the external potential. 

■ ' 1. Introduction 

In this paper, we study the colhsion of two fast sohtons in the presence of a 
(time-dependent) external potential that varies slowly in space compared to the 
size of the solitons. We show, for a class of typical local and nonlocal nonlinear- 
O . ities, that if the initial relative speed of the solitons is \\v\\ ^ 1 and the spatial 

variation of the external potential is sufficiently slow, then the solitons pass each 
t:;;;}^ . othcr aluiost blindly: The L^-norm of the difference between the true solution 

O I and the one corresponding to a configuration of two solitons moving in the ex- 

(^ ! ternal potential decays algebraically with ||w||, up to times of order log \\v\\, after 

the collision. This is an example where the solitary waves for NLS display both 
their "wave" and "particle" nature. They pass each other almost blindly because 
Jh I they are localized waves with high relative speed and relative phase, while their 

center of mass dynamics is approximately that of a classical particle in a spatially 
adiabatic external potential. 

The problem of asymptotic behaviour of multi-soliton configurations (scatter- 
ing theory) for the nonlinear Schrodinger equation without an external potential 
has been addressed in [1] and [2]; see also [3]. In these papers, the authors 
prove, under rather strong spectral assumptions on the linearized equation, the 
asymptotic stability of multi-solitons in three (or higher) dimensions. The main 
ingredient of their analysis is asymptotic stability of single solitons and dispersive 
estimates (which are related to the "charge-transfer model"). Here, our results 

and approach are different: We study the long-time dynamics of the collision 
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of fast solitons in the presence of an external potential rather than the asymp- 
totics, and we use softer yet more robust techniques that allow for treating a wide 
class of systems under weak assumptions. Furthermore, our analysis holds in any 
dimension A^ > 1. 

There has been considerable progress in understanding the long-time dynamics 
of single solitons in spatially adiabatic external potentials and in the presence of 
nonlinear perturbations, [H El [6l [71 [H [9] . The analysis below together with addi- 
tional mild spectral assumptions can be extended to study the effective dynamics 
of multiple solitons with low velocities in slowly varying external potentials (and 
in the presence of nonlinear perturbations) as long as the soliton centers of mass 
are well separated. \j 

1.1. Description of the problem. We consider the nonlinear Schrodinger equa- 
tion 

(1) zdMx, t) = (-A + \4(x, t))ij{x, t) - f{ij{x, t)), 

where A = X]i=i ^ i^ ^^^ A^-dimensional Laplacian, with A^ > 1, V^ denotes 
the (time-dependent) external potential, with 

Vh{x,t) = V{hx,t), 

and / is a focusing nonlinearity 



such that /(^) = fiifj). 

We now discuss the various assumptions we make, which are simultaneously 
satisfied by typical local and Hartree nonlinearities, see Remark [1] below. 

(Al) Global well-posedness. The nonlinear Schrodinger equation dTj) is globally 
well-posed in H^. 

We refer the reader to [13], chapter 6, for well-posedness of ([1]) in energy space 
for time-independent potentials, and [7] for the case of time-dependent external 
potentials and nonlinearities. We make the following assumption on the regularity 
and symmetries of the nonlinearity. 

(A2) Nonlinearity. Let F : H^ ^ M be the functional such that its Frechet 
derivative F' = f. We assume that F G C'-^iH^; R) and that F(T-) = F(-), 



We note that for the generahzed KdV equation, there has been some recent progress in 
understanding the coUision of a fast thin sohton with a slow broad soliton in the absence of an 
external potential, see [inilll]; and also [T^] for a recent review about problems related to the 
stability of solitons. 



where T is a translation 

a rotation 

T^ : u{x) -^ u{R^^x), R G SO{N), 
a gauge transformation 

T^ -.uix) ^e'^u{x), 7G [0,27r), 

or a boost 

We are interested in the dynamics of multi-sohtons, so we assume the existence 
of sohtary wave solutions when V = 0; see for example [I3], chapter 8, for a 
discussion of solitary waves for NLS. 

(A3) Solitary waves. When V = 0, there exists an interval / C M such that, 
for all yU G /, ([1]) admits solitary wave solutions of the form 

where 

a = (a, V, 7, ;u) G M^ X R^ X [0, 27r) x /. 

Here, r/^ is a positive and spherically symmetric function satisfying the 
nonlinear eigenvalue problem 

(2) (-A + ^)7y,-/(r/^)=0, 

(3) |||xP?7^||l2 + |||xnV?7^|||L2 + |||xp9^?7^||L2 < oo,VyU G /, 



and 



Let 



?7^ oc e ^ll^ll as ||x|| -^ oo. 



m(yu) = - I dx rf 



2 ' """" ■'^' 
the "charge" of the soliton. We assume that 

d^m{fi) > 0, 

which implies orbital stability, see p! ^fT5l flG] . 

We require some local properties of the nonlinearity, which are satisfied for 
classes of local and nonlocal nonlinearities, see Remark [H 
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(A4) Localization. We assume that 

and 

where 77^^ = M(o,o,o,a«i)' ^ = 1, a; or 9^.. Here, ^ G (0, min(y/ir, ^/I^)) and 



C > are constants that are independent of 02 and V2. 
(A5) For g G Lp(]R^), p > 1, m = X]"=iMo-i, where u^ appears in (A3), and 
w G H^ with 11^11^2 < 1, we have 

\{gJ{u + w)-f{u)-f{u)w)\<C\\w\\l,, 

where C is a constant the depends on g and fii,i = 1,- ■ ■ ,n. 

We make the following assumption on the external potential, which, among 
other things, guarantees well-posedness of ([T]) in H^, inspite of the fact that the 
energy in no more conserved, see [7j. 

(A6) The external potential V G W^''^(R; C^{R^)). 

We now discuss the initial condition. We are interested in the collision of solitons 
with high relative speed. A 2-soliton configuration plus a fluctuation is given by 

(4) ij{t = 0) = (p{x) = e5^^-%-j^(x -Si) + ei^^"'r]ji^{x-a2) + w, 

with 0,1,0,2, Vi,V2 G M^, w G H^, and Jli,Jl2 G Jq, where Jq C I\dl is a bounded 
interval such that its closure Jq C I\dl. We assume that 

11"^! — "^211 ^ ( inf m\jj))~^ 

with pi^ = 0(1). 

We assume that the fluctuation w is small. More specifically, w G H^ such 
that 

||U'|||2 < C\\Vi — ■^211""^- 

We did not impose any condition on the directions of the relative speed and 
position of the solitons. In particular, we can have 

(ai - 02) ■ {vi - V2) < 0, 

which is the case corresponding to colliding solitons. We remark later how one 
obtains better estimates in case the solitons are escaping each other. 

In what follows, we denote by fo := ^1 ~ ^2, the initial relative velocity of the 
solitons. 



1.2. Main result. We are in a position to state our main result, whose general- 
ization for fast n-solitons, n > 2, is straight forward. 

Theorem 1. Consider the nonlinear Schrodinger equation (J\) with initial condi- 
tion given by Oj, and suppose assumptions (A1)-(A6) hold. Then, for any fixed 
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a G (0,1), ||fo|| ^ (iiifM6-fo "^'(/^))~^^^ ^'^^ ^ ^ (ii^fMS/o "^'(/^))^^^; ^^^ solution 
of the initial value problem can be written as 

iP{x,t) = e*^i*+^^-(^-"^)r/^,(x - ai) + e'^^'+^''^<^-''^^r]^^{x - as) + w{x,t), 

for all tE [0, Tq,), Tq := Camin(log llfoll, 2| log/i|), and 

(5) sup \\w\\L2<C'i\\vor'^ + h'~n, 

te[o,T„) 

where the constants C,C' > are independent ofvo,h and a. Furthermore, the 
parameters a^, Vi, 7,, /ij, i = 1,2, satisfy, for t G [0, r^), the following equations 

dtai = Vi + O(||t;o|r('~") + /i'('-") + e-«ll"^-"2"), 

dm = -2VVh{a„ t) + O(||i;oir^'"°^ + /^'^'""^ + e-^"'^-'^^"), 

dtl^ = /i. + f - H,(a„ t) + 0(|K;o|r(^-") + /i'^^"") + e-^""-"^!'), 
/or some ^ G (0,niin(y^, JJ12)) that is independent of \\Vq\\ and h. 



In particular, for ||uo|| ^ 1 sufficiently large, and h = 0(||t;o||~^)) the solitons 
preserve their shape, in L^-norm, up to times log ||vo|| after the collision, such that 
the dynamics of the centers of mass of the solitons is approximately determined 
by the Hamilton equations for two classical particles in the external potential. 

Our analysis relies on three main ingredients. First, using a skew-orthogonal 
(or Lyapunov-Schmidt) decomposition property (Proposition [H Sect. H]), we de- 
compose the solution of ([1]) with initial condition close to a 2-soliton configu- 
ration, as described by (jl]), into a path belonging to a symplectic manifold of 
2-soliton states, and a part describing a fiuctuation skew-orthogonal to the mani- 
fold. The dynamics on the 2-soliton manifold is obtained by the skew-orthogonal 
projection of the Hamiltonian fiow generated by the nonlinear Schrodinger equa- 
tion in a small tubular neighbourhood of the 2-soliton manifold onto the latter 
(Proposition [2], Sect. [5]). As for the fluctuation, we control its L^-norm using 
charge conservation and skew-orthogonal decomposition (Proposition [3l Sect. [6]). 
The main difference between our approach and the one for studying the effective 
dynamics of a single soliton in an external potential, as for example in [4], is 
that we control the L^-norm of the fluctuation using charge conservation, rather 



than controlling its iJ^-norm by using an approximate Lyapunov functional and 
proving constraint positivity of the Hessian, Eq. ( fT3l) below, under additional 
assumptions that are verified in the case of special local nonlinearities. Unlike 
the L^-norm, the H^-xioim of if) grows like ||fo||) s-nd we lose control over Hu^Hi^i 
as ||yo|| -^ oo. 

Remark 1. We now give some concrete examples for which assumptions (Al)- 
(A5) are simultaneously satisfied. 

An example where assumptions (Al) - (A3) are satisfied is when f is a Hartree 
nonlinearity , 

such that W is positive, spherically symmetric, belongs to L^ + L°°, with p > 
Y,p > 1, and decays at infinity, W —>■ as \\x\\ —>■ cxo; see [I3l[7]. The localization 
property, assumption (A4), is satisfied if in addition W decays exponentially fast. 
We now verify that (A5) holds for p > 2. From the form of the nonlinearity, we 
have 

\{g, f{u + w)- f{u) - f'{u)w)\ < C{\{g, {W^\w\^)u)\ + \{g, {W ^ \uw\)w)\) . 

Applying Holder's and Young's inequalities, we have 

\{g,{W^\w\'')u)\< WguW^AW^H^L'^ 

< ii5'iiL'j'ikiu°°iii^iU''iiknui 

n 

< llfi'lL'j'lX] lkilU°°)II^IU''lklli2 

<C\\w\\l,, 
where q = p or oo (W G L^ + L°°) and 1 = 1/q + 1/q'. Similarly, 



{g, {W -k uw )w) < gW * 


Uw\\\l2\\w\\l2 




\\W -k {\u\\w\) LI W L2 


< \9\ ^ \ 


\W \li \uw\ lA ^ \l2 


< 1 ^ , % 


\W\li\\u\l2 W\l2 


< \g\ 2q 


n 
i=l 


<C\\wfr2 





Therefore, assumption (A5) is satisfied. 



Another example where the various assumptions are satisfied is when f is a 
local nonlinearity. For example, (Al) and (A2) are satisfied if f is of the form 

where he C^iM+.M) with 

d^h{r) < C{1 + r"^^'), ^ = 0, 1, 2, 

a E (0,;^^32)' ^ ^ 3, and a G (0, cxd) if N = 1,2; see for example [131 El for 
a discussion of well-posedness in H^. Solitary wave solutions appearing in (A3) 
exist, if, in addition, 

— CX3 < lim h{r) < /i 

r— >oo 

— oo < lim r~"^h{r) < C, 

r— >oo 

and there exists ro > with 

ro 

h{r)dr > fiVQ, 



see [T71IIS]- The condition of orbital stability can to be checked for each nonlin- 
earity, see [HI ESI [IS] ■ Assumption (A4) follows directly from (A3) and the form 
of the local nonlinearity. Furthermore, assumption (A5) is satisfied if 

'2k — 1 J 

sup r^^dj.h{r) < oo, k = 1,2. 

An explicit example of a local nonlinearity that satisfies all the above hypotheses 
is 

where Xe,s, 6 ^ I, is a smooth regularization which is chosen such that (A5) is 
satisfied. For example, 

Jl, if \y\'9<'-^) < e/2 

\\y\'-\ if |y|^^"(^-i' > ^ ■ 

More generally, f can be a sum of both local and nonlocal nonlinearities. 

Remark 2. We now remark on special cases where one can obtain a control of the 
fluctuation over different (and longer) time scales. Assume (Al)-(A6) hold, and 
suppose, for the sake of simplicity, that h = 0, which corresponds to a spatially 
flat potential. 

(1) Large separation. If the soliton centers of mass are initially separated by 
a distance d ^ max(^L, ^L, | loginf^g/,j m'(yu)|) and H'^iH, ||t?^2|| = 0(1), 



with \\w\\l2 = 0{e~^'^) for some x > 0, then one obtains a result similar 



to TheoremUl such that sup^grQ,^^) ||w||i2 < C/y/d, for any e G (0, 1) and, 
forte[0,d^), 

a,{t)=a, + tu, + 0{d-^'~'^) 
u,it)=v, + Oid~^'-'^) 
fi,{t)=Jl, + 0{d-^'"'^) 

7.(t) = 7. + tiJi, + J - V^(0, t)) + 0(rf-(i-^)), 

with i = 1,2, see Sect. 0. 
(2) Escaping solitons. Suppose that the solitons escape each other with a high 
relative speed 

(ai -02) ■ {vi -V2) > 0, 

and \\vq\\ ^ (inf^g/g m'(/i))~^, ||w^||l2 = 0(e~^l'^'"l') for some x > 0, then, 
for any fixed e G (0, 1), there exists a contant C, independent of \\vq\\ and 
e, such that 

sup \\w{t)\\L2 < C\\vo\\-'^ 

and, for t G [0, H^oin, 

ai{t)=Zi + tVi + 0{\\v4-^+') 
v,{t)=v, + 0{\\v,r^+^) 
/i,(t) = /I, + 0(||t;oir'+^) 

7.(t) = 7. + t{J^^ + J - V{Q, t)) + 0{\\v,r^+^), 

with i = 1,2, see Sect. 0. 

The organization of this paper is as follows. In Sections [2|, we recall some basic 
properties of the nonlinear Schrodinger equation. In Section [3l we recall the 
soliton manifold, and we introduce the 2-soliton (or, more generally, ra-soliton) 
manifold. In Section |U we prove the skew-orthogonal decomposition property 
for elements of neighbourhoods in H^ that are close in (L^-norm) to a two- 
soliton manifold, which is a central tool in our analysis. In Section [5], we use 
the skew-orthogonal property and the nonlinear Schrodinger equation ([T]) to find 
the reparametrized equations of motion corresponding to the parameters on the 
two-soliton manifold, and in Section [6l we control the L^-norm of the fluctuation 
using charge conservation and the skew-orthogonal decomposition. In Section [7] 
we prove Theorem [1] by combining the results of Propositions [H |2] and [3l We 
finally remark on separating solitons in Section [HI 



1.3. Notation. 

• In the following, L^i^I) denotes the standard Lebesgue space, 1 < p < oo, 
with norm 

\Lo^=esssup{\f\), feL^il). 
We denote by (■, ■) the scalar product in L^(]R^), 

{u,v)=Re uv, m,wgL2(M^). 

Given the multi-index a = (ai, ■ ■ ■ , un) G N^, we denote |a| = J2i=i ^i- 

Furthermore, 9" := 9°^ ■ ■ ■ 9°^. 

For 1 < p < oo and s G N, the (complex) Sobolev space is given by 

PF^'P(M^) := {u e S'{R^), dy e LP(M^), |a| < s}, 

where iS'(R^) is the space of tempered distributions. We equip PF^'*' with 
the norm 

\\u\\ws,v = ^ WdyWip, 

a,\a\<s 

which makes it a Banach space. We use the shorthand W^'"^ = H^. 
Given / and g real functions on M^, we denote their convolution by *, 

f-kg{x) := / dy f{y)g{x-y). 



1.4. Acknowledgements. W.A.S. thanks Catherine Sulem for pointing out ref- 
erences [IOl[ni[T2]. 

2. Hamiltonian structure of the nonlinear Schrodinger equation 

In this section, we recall some basic properties of the nonlinear Schrodinger 
equation ([T]), see for example [16111]. We will use these properties in the following 
sections. 

The space i/^(M^,C) has a real inner product (Riemannian metric) 



(6) (u, v) := Re / dx uv 

for u,v E H^iM^ , C). □ It is equipped with a symplectic "form'' 

(7) Lu{u, v) := Im / dx uv = {u, iv). 



^The tangent space a.i ij) e H^ is T^,H^ = H^. 
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The Hamiltonian functional corresponding to the nonhnear Schrodinger equation 
(HD is 

(8) Hy{tP) := i I \Vij\'dx +^-lv\tlj\'- F(^). 

Using the correspondence 

H\R^, C) < — > H\R^, R) © H\R^, R) 

i^^ < > J, 

where J := j n ) ^^ ^^^ complex structure on if^(M^,]R^), the nonlinear 

Schrodinger equation can be written as 

Furthermore, 



(m, v) = dx (Reu liau) ( ^ ) 
J \lmvj 



Rev 
Imw 



We note that since the Hamiltonian functional Hy defined in ([8]) is nonau- 
tonomous, the energy is not conserved. For tp G H^ satisfying ([1]), 

dtHv{^) = ^JdxidtVM\ 

see [7] for a proof of this statement. Still, Hy is invariant under global gauge 
transformations, 

and the associated conserved quantity is the "charge" 

N{^):=^jdx\i:\'. 

The assumption d^m{jji) > implies that tj^ appearing in assumption (A3) is a 
local minimizer of Hv=o{'ip) restricted to the balls Bm := {^ G H^ : N{tp) = m}, 
for m > 0; see [T^ [15] . They are critical points of the functional 

(9) £,{^P) ■■=\j dx (I VV'P + ^i\i^?) - F(^), 

where /i = ^{m) is a Lagrange multiplier. 
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3. SoLiTON Manifolds 

In this section, we recall the definition and properties of a single soliton mani- 
fold (see in El [6l [71 [HI [9]), and we introduce the multi-soliton manifold. 

3.1. Soliton Manifold. We introduce the combined transformation Tav-y, which 
is given by 

i^av-y := Ta^yiJ = e'^"^^<^-''^^'^H{x - a), 

where w, a G M^ and 7 G [0, 27r). We define the soliton manifold as 

Ms := {r/a := T,„^r/^, a = (a, v, 7, ^i) eM."" xR"" x [0, 27r) x /}, 

where / appears in assumption (A3). If /'(O) = 0, where / appears in ([1]), then 
JCM+. 

The tangent space to the soliton manifold TW^ at 77^ G TW^ is given by 

%,^Ms = span{Et, Eg, Eb, Es}, 

where 

Et := VaTf r/J,=o = -Vr/^ 



In the following, we denote by 



E, 


= 


d.T^V, ,- 


=0 — 


iVi^ 


Eh 


= 


2V.T,^^ 


v=Q 


= iXTj^ 


Es 


= 


df^Vi^- 






8 by 








Cj := 


-dx, , j = 


= 1,- 


■■,N, 


Gj+N 


■ iXj, J - 


= 1, 


■■■ ,N 


e2N 


+1 


■=h 







(10) e2N+2 ■■= d, 



M' 



which, when acting on rj^ G A4s, generate the basis vectors {eaVo-}a=i'^ oi%i^M.s- 
The soliton manifold A^s inherits a symplectic structure from {H^,uj). For 
a = (a, V, 7, /i) G M^ X R^ X [0, 27r) x J, 



n^ := P^J-'P. 



-1 ^ 



where Pa is the L^-orthogonal projection onto T^^TW^. 

We have the following easy lemma, which we prove in the Appendix. 

Lemma 1. If d^m{^) > 0, then Q„ is invertible. 
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Explicitly, we have 






^cr %j„Ms '■ = {(ea?7^,'ie;3?7^)}i<Q,/3<2Ar+2 






(11) 


/ -m{lJ,)lNxN 

m{n)lNxN 








-lfm'(/i)\ 
am'{fi) 

/ 



where Inxn is the N x N identity matrix, and {■)'^ stands for the transpose of a 
vector in M^; see the proof of Lemma [1] in the Appendix. 

3.2. Group structure. The anti-self adjoint operators {ea}a=i,---,2N+i defined 
in fITOl) form the generators of the Lie algebra g corresponding to the Heisenberg 
group H^^^^, where the latter is given by 

{a,v,j)-{a',v',i) = {a",v",j"), 

with a" = a + a' , v" = V + v' , and 7" = 7' + 7 + |w ■ a'o Elements of g satisfy the 
commutation relations 

(12) [cj, Cj+at] = -e2N+iSij, i,j = !,■■■ ,N, 
and the rest of the commutators are zero. 

3.3. Zero modes. The solitary wave solutions transform covariantly under trans- 
lations and gauge transformations, i.e., 

for all a G M and 7 G [0, 27r). There are zero modes of the Hessian, 

(13) C,:=-A + fi-f'{r],), 
associated to these symmetries. We have the following lemma. 
Lemma 2. 

with {iC^YX = for any vector X G T^^TW^. 

Proof. Differentiating £'^{T^ri^) = with respect to a and setting a to zero 
gives 

(14) S"{r]^)VaV,.{x - a)\a=o = C^Et = 0. 

Similarly, differentiating £' {T^rj^) = with respect to 7 and setting 7 to zero 
gives 

(15) C^Eg = 0. 



This structure was noted for the case TV = 1 in [5]. 
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Using ([2]), we have 

(16) C^Eb = (-A + /i - f'{T]^))ixr]^ = -iV^r]^ = iEt, 
Furthermore, differentiating ([2]) with respect to fi gives 

{-A + i2-f'{r],,))Es + v>. = 0, 
and hence 

(17) C^E, = i{ir]^) = iEg.D 

3.4. Two-soliton manifold. We now discuss the manifold corresponding to two 
sohtons. It is given by 

Ml := {(Va^Va^), (^^ = {ai,Vi,-f^,fi^) G M^ X M^ X [0, 27r) xl, i = 1,2}. 

The tangent space to Ai^ is 

nv.„,.,)M! = {(Xi,X2), X, G T,^Ms, t = 1,2}. 

We introduce the embedding mapping 

E : A?2 ^ H\ 

whose action on A1^ and TAi^ is given, respectively, by 

E(^ai,^a2) = r]a^ +Va2 e H\ 

E(Xi, X2) = Xi + X2 G Ti/^ + TiJ^ 

In what follows, M'^ and TM^ denote E{Ml) and E(TA^^) respectively. 
4. Skew-orthogonal decomposition 
Let I be the same as in assumption (A3). We define 

S := {a = (a,w,7,/i) G M^ x M^ x [0, 27r) x /}, 
and let 

E° := {a = (a, w, 7, /i) G M^ x M^ x [0, 27r) x Jq, with Jq C /\9/ bounded}. 
We define 

^d,K •= {('^15 ^2) G S° X S°, ||ai — 02!! > (i or \\vi — V2\\ > k,}- 

In other words, for (cri,a2) G S^^, the centers of mass of r/g-i and r/o-j are either 
separated by a distance larger than d or their relative speed is larger than k. 
We consider the neighbourhood Us^d,K C H^ defined by 

U5,d,K ■= {ip e H^, sup Wijj - 7]^^ - ?7<x2||l2 < ^}. 
We have the following proposition. 
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Proposition 1. Suppose (A2) and (A3) hold. Then, ford <^ inf^g/^ m'(/i) and 
^ ^ inf^s/o"^'(M) ('^^ ^ ^ "^^^'^^if ' ;if ' I logi^f^e/o "^'(/^)l)^. ^^ere exist unique 



sttc/i i/iai 
(18) 
and 
(19) 

/or a// Xj G 7;,^^(^,)A^^, 2 = 1,2. 
Proof. We define the mapping 



■0 = Vaiii,) +Va2ii>) +W, 



uj{w,Xi) = 0, z = 1, 2, 



G : t/5,d,K X S^ 



D4Ar+4 



by 



G'(z^,(cri,cr2)) : = 



/ uj{i> - r]„^ - ?7^2, ei?7ai) \ 






\w(?/; - r/^j - r/^2, e2N+2V<T2)/ 
Tlien (TT8l) and (TTOll are equivalent to (cti (■?/;), a"2('?/')) satisfying, for a given ■?/; G 
Us^d,K, the equation 

G(V^,(ai(^),a2(V^))) = 0. 

We use the implicit function theorem to show that there exist unique o"i(^), 0'2{ip) G 
E such that G{ip, {ai{ip), <J2{ip))) = 0. 
First, note that, by construction, 

(20) G(r/., + r/.„(ai,(T2)) = 0. 
Furthermore, 

(21) GeC\Us,,,.xj:l,;R'^+'), 

since it is hnear in i(j and rjcr^, i = 1, 2, and it is differentiable in ai G S", z = 1, 2. 
We still need to show that 9(o.j_o.2)^('?o-i +Vt2j (3^1) ^2)) l^i =0-1,52=0-2 is invertible for 
K, 3> ^-r — ^—n-T (or d S> max(^L, ^=, ^-f — ^ ,, , )). 
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We know that the matrix 

{uj{eaVa,ef3Va)}l^f3t\, 

is invertible, see ( II ip . Lemma [H 
We write 

(22) e^r]^,J^^, =: e^^'-^'-^'^^h^fsix), «, /? = 1, ■ ■ ■ , 2iV + 2, 

which corresponds to a decomposition where the fast oscillating term (in space) 
gf {vi--«2)-a; jg separated from the slowly oscillating term (in space) hap- Let \\vm\\ '■= 
max(||wi||, 11^211, 1). It follows from the fact that f E C"^ (assumption (A2)) and 
the exponential localization in space of the solitons (assumption (A3)), that there 
exists ^ e (0, min(y^, y^)), which is independent of \\vi — V2\\, and a constant 
C that dependends only on yUi and ^2^ such that 

(23) II Vlk3.i(iRiV) < C|K;„f e"«ll"^-"^ll, 

for a, /3 = 1, ■ ■ ■ , 2A^ + 2o Suppose that k^ 1. Let v := Vi — V2. Using that 

where 

V 

L := -2i^ ■ V,, 

IK'll 

and integrating by parts three times, we obtain 

t^(ea?7ai, 6/377^2) = / {L^e^'"''')K(i{x) dx 



(24) = / e^-"(L*)3/i„^(x) dx. 

Moreover, 

(25) \\{L*fhai3\\L^ < ||t;|r^||/i„/3||H/3,i. 

Eqs. ([221) - (ESD yield 

\uj{eariai,epr]^^)\ < C||t;f ||t;||~^ 
<C||t;||^^ 

(26) < Ck-\ 



%ore generally, if / G C''{H^,H-^), 

l|/ia/3||w-+i.MK")<C|lz;™fe-«ll'^^-"=". 

For example, in the case of local nonlinearities, the above estimate holds for any r > 1, in which 
case we obtain better estimates. 
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(Suppose alternatively that d ^ max(^L, ^L) with ||fo|| = 0(1) fixed. Then 
it follows from fl23l) that 

|c^(ea^<xi, 6/377^2) I < Ce'^"^, 

for some positive constant C that depends on /ii and ^2 and ^ G (0, min(^y^, ^JJI^)). 
Hence, for 



-') 



and 



5 <^ inf m'^jj) 



K :» ( inf m'(n)) ^ 



(or (i > max(--|=, --^, | loginf^g/o m'(/i)|)), the {AN + 4) x (4iV + 4) matrix 

is invert ible. 

Invertibility of 9(o-i,cr2)G', together with fl20|) . fl2Tl) and the implicit function 
theorem, g imply that there exist unique C^ maps cti (■?/') and cr2(^) such that 

G(^,(ai(^),a2(^))) = 0. 

D 

Remark 3. The group element Tav-y G H^^^^ is given by 

rp _ p-a-d^ i^ i-y 

/t follows from fTB) that T~^YTav-y E g if Y E g. Furthermore, it follows from 
translational invariance that uj{TavyU,Tavyv) = w{u,v), for u,v G L"^. Therefore, 
we have from Proposition Ul that 

uj{w, Y{ri^, + r/<,J) = uj{w', Y'{r]^'^ + r/^^)) = 0, 

VF G g, where Y' = T'J^YTav-y e g, w' = T^^J^w, and r]^' = T^JyVa- 

5. Reparametrized equations of motion 

In this section, we apply the skew-orthogonal property to obtain reparametrized 
equations of motion for the parameters that characterize the projection of the true 
solution of ([T]) with initial condition onto Ai^. 

We assume that the hypotheses for the skew-orthogonal decomposition. Sect. 
m hold. We will verify in the proof of the main theorem that for large enough 
lluoll and small h, this is indeed the case over a certain time interval. 



See for example 
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Proposition 2. Consider ([7]j with initial condition Q), and suppose that (Al)- 
(A6) hold. Assume further that there exists r > such that, for t G [0,r), 
^/^(t), the solution of (QP with initial condition (p, is in Us^d,K, where 6 is given 
in PropositionUl Then, for \\vq\\ ^ 1, there exists a positive constant Cq < 1 
independent of \\vo\\ and h, such that, for ||u^||l2 < Co, the parameters a, = 
{ai,Vi,'~fi, fii), i = 1,2, satisfy the equations 

(27) dta, = v, + 0{\\w\\l, + h^ + e-«ll«i-'^^ll), 

(28) dtv, = -2VaMai, t) + 0{\\w\\l2 + h^ + g-^H"^-"^'!) 

(29) ^a^ = f^^ + J- yh{a^, t) + 0{\\w\\l, + h' + e"«ll'^-«^ll), 

(30) dti^i = 0{\\w\\l2 + h'' + e-^ll'^i-'^^li)^ 

for some constant ^ G (0,min(y^, a//U2)) that is independent of \\vq\\ and h. 

In what follows, we denote by C a positive constant that is independent of \\v\\ 
and h, but that may change from one line to another. 
Proof. We first find the equation of motion for 

Using Proposition [H we have 

(31) ui = r]f,^ + 7]^'^ + w' , 

where r^^/ = T-^l^^^7]„^ and w' = T~^l,^^w. Here, a'^ = Oa - Oi, v'^ = V2 - Vi, 
72 =72 — 71, and /i'a = /i2- It follows from Remark [3] that 

(32) uj{w',Xi+X2) = 0, 

for all Xi G %^^Ms and X2 G T;,^, A^s- 
We define the coefficients 

Cj := dtaij - vij, j = !,■■■ ,N, 

CN+j ■■= -^dtVij - Va,Vh{ai,t), j = !,■■■ ,N, 

C2N+1 ■= /ii - jvl + -dtai ■ vi - Vhiai, t) - %, 

(33) C2Ar+2 := -<9t/i. 

Note that 

.,2 

(34) 6-5(^1 -^+271 )A7/^(x + ai) = Ami + ivi ■ Vmi - ^ 

(35) e-i('^-^+270/(^(^ + a,)) = f{u,). 
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Differentiating ui witli respect to t and using ([1]), fl33l) - fl35l) . we get 



2N+1 



(36) dtui = -i((-A + fii)ui - f{ui)) + ^ c„e„ni - iTZyUi, 

Q!=l 

wliere 

■^v = Vh{x + ai, t) - Vh{ai,t) - x ■ Wh{ai, t). 
In otlier words, 

2N+1 

(37) 9iMi = -i£:^^(Mi) + ^ c^eaMi - ilZvUi, 
wliere £^ is defined in (^. Recall tliat 

wliicli implies 

(38) S'^^im) = C^,{r]^> + w') + N^.iv., + w'), 

wliere 

>Cmi = (-A + /ii - /'(r/^J) = ^^',(^/.J 
and 

iVMi(^4 + ^') = /(^w + ^4 + ^') - /(^Mi) - /'(^m)(^4 + ^')- 
Substituting ( 1311) and ( l38l) into ( 1371) . we obtain 

2N+1 2N+2 



dtw' ={-iC^^ + ^ c„e„ - inv)w' + N^,^{ri„i^ + w') + ^ CQe„?7^^ - ^7^v'^A'l 

o=l o=l 

27V+1 

(39) - (9(77^^ + (-i^^i + ^ c„ea - iTZv)r]„'^. 



a=l 



To obtain the equations of motion for ai, f i, 71 and /ii, we use the skew-orthogonal 
property to project (1391) onto 7^^ TW^- 

It follows from ([32]) that (iw',X) = for all X e T^^^Ms- Therefore, 

(40) dt{tw\X) = dtfii{zw',d^,X) + {idtw\X) = 0. 
Substituting the expression for dfw' given by fl39|) in fHOl) . and using 

(41) e: = -e„, « = l,---,2iV + 2, 
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we have 

2N+2 2N+2 

Q=l «=1 

(42) 

2N+1 
o=l 

Some of the terms in the above equation drop-out due to the zero modes of the 
Hessian. It follows from (fT4l) - (TT7j) . Lemma [2], that 

X = i^fi^X e Trj^^Ms if -^ e Tri^^Ms: 

and hence 

{C^y,X) = {w',C^,X) = -uj{w,X')=0. 

Together with (gl]) and (iSl), this yields 

2N+2 2N+2 



y^ CaUj{ear]^^,X) =(7^y?7^,,X) + ^ Cq {ieaw'.X) + {nvw',X) + {iN^^{'q^>^ +w'),X) 

Q=l a=l 

2Ar+l 

(43) + (7^yr7,^ , X) + (zSiT/,. , X) + (r/,. , [C^, + ^ ^ c,e„)X) . 

a=l 

We now estimate each term appearing in the right-hand-side of (H3l) with X = 
epTj^^, /3 = 1, ■ ■ ■ , 2X -|- 2. Note that it follows from assumptions (A3) and (A6) 
that 

Wnve^Tl^, IU2 = 0(/i2), /3 = 1, ■ ■ ■ , 2X + 2, 
and from (A3) that 

II^IIl^ = l|e/3r/Mi lU^ = 0(1), /3 = 1, ■ ■ ■ , 2X + 2 
||e„X||i2 = lle^e^r/^, ||l2 = 0(1), a, /3 = 1, ■ ■ ■ , 2X + 2. 

Hence, Holder's inequality, (A3), (A6) and the fact that V is real yield the esti- 
mates 

(44) \{nvV,.,X)\ = Kr/^„7^yX)| < h.AWWT^vepV.AW < Ch' 

(45) \{nvw',X)\ = \{w\nvX)\ < \\nvef,r]^,\\L2\\w'\\L2 < Ch^\\w'\\L2 

(46) \{nvVa',X)\ = |(r/^.,7^yX)| < \\T]^dL4'^ve(3V^.^\\L^ < Ch\ 
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We also have from (A3) and Holder's inequality that 

2N+2 2N+2 

I ^ c„(ie„w',X)| = I ^ c^{iw\eaX)\ 
(47) <C||c||||w;'||i2 

where ||c|| := maxQ,=i^... ^2Ar+2 \ca\- 

We now use assumptions (A3)-(A4) to evaluate \{iN ^-^{tj^i^ +w'),X)\. It follows 
from (A3) that iX = iepT]^^ G W , P > 1, which, together with (A5), yield 



l(^(/(^m + Va' + w') - f\ri^, + r]^>) - f\ri^, + iq„:)w'),X)\ < C\\w^ 



'\\2 

Il2- 



It follows from the boundedness and the exponential localization of the solitons 
in space, (A3), and the fact that / e C^, (A2), that 

||H/'(r/,Jr/.,,X)| < ||/'(r/,JlUoo||ry,,X|Ui 

for some ^ G (0, min(y//i7, y^)) which is independent of \\vq\\ and h. Moreover, 
it follows directly from (A4) that 

and 

mf'{v,^+v.O - f'{v,J)w',x)\ < \\w'h4if'iv,.+v.O - f'iv,J)x\\L^ 

Therefore, 

(48) \{zN,,{rj^,+w'),X)\ < C{\\w'\\l^ + e-^^^^^~'^-^^). 

To evaluate the remaining terms, we use the fact that 77^^ and rj^i are expo- 
nentially localized in space, while their relative fast oscillating phase is 

||t;|| = llt;^ -17211 > |K'm,||, 

where \\vm\\ := max(l, ||t;i||, ||t;2||). 

When estimating an upper bound for |(«9i?7o-^, 6/377^^) |, the partial derivative 
with time contributes ||fm|P, since 



(49) dtVa 



N , 2N 



y^ dtajCj + 7; ^ (^tVjCj + {da + ——)e2N+i + dt^e2N+2 



2 

j=l j=N+l 



Va- 
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However, using fl22l) and fl23l) . and integrating by parts twice in space, we can pull 
a factor of \\vi — f2||~^ from the fast oscillating term e^^^'^-^'O-a;^ ggg ^^j^g discussion 
below (l23l) in the proof of Proposition [H Hence 

(50) mv.',,epv,.)\<Ce-^^\^^~^-K 
Furthermore, (A2) and (A3) yield 

(51) |(ry..,£,,e^r/,,)|<Ce-«ll---ll. 

Again, using (l22l) and fl23|) and integrating by parts twice in space to pull a factor 
of ||fi — f2||~^ from the fast oscillating factor e^^^^~'"^^'^ , we have 

2N+1 

(52) |(r/4,(2 Yl Caeje^r/^,)! < C\\c\\\\v\\~'e-^^^-^-'^'K 
From (03]) - ([52]), we have 

2N+2 

(53) I ^ c„c.^(e„r7^,,e^r/^J| <C[|K/7||i. + ||c||(||«7|U2+|K;||-2) + /,2 + e-«K-a2ll]^ 

a=l 

for /5 = 1, ■ ■ ■ , 2A^ + 2, where we used 

due to translational invariance. 

Using Lemmadl (ITT!) and (1531) . and assuming ||w||l2 and ||t'||^^ < ^||r2^J|, we 
obtain the estimate 

||c|| <C(||u;||i2 + /i2 + e-^ll"i-"2l'). 

Recalling now the definition of Cq, a = 1, ■ ■ ■ , 2A^ + 2 (see fl33l) ). we conclude fl27j) 
- (EO]), with i = 1. 

To get the equations of motion for 02, f2, 72 and /i2, we consider U2 = T~J^^^^ip, 
and we repeat the above analysis with 1 <-> 2. D 

6. Control of the fluctuation 

We now control the L^-norm of the fluctuation w using conservation of charge, 
the skew-orthogonal property. Sect. IH and the reparametrized equations of mo- 
tion. Sect. [51 

Proposition 3. Consider (J\) with initial condition Q), and suppose that (Al)- 
(A6) hold. Assume further that there exists r > such that, for t G [0,r), 
il){t) G Us^d,K, where 6 is given in PropositionUl Then, for \\vq\\ ^ 1 and h <^1, 

sup lk(t)||i2<c"(|K;oir^+" + /i'(^"")), 

te [0,Ca min(log ||i)o || ,2| log /i|)) 
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for some positive constants C and C that are independent o/vq, h, and a G (0, 1). 
Proof. From conservation of charge (L^-norm) of the solution of ([T]), 

and skew-orthogonal decomposition (Proposition [T]), we have 

(54) ll^lli. = ll^lli^ + \\ri^,\\l, + \\ri^,\\l, + 2Re{r]^,,r]^,) = H^Hi^, 

where we used 

{w,Vaj) = -cj{w,iri^^) = 0, 
and 

for j = 1,2. 

Differentiating (!54l) with respect to t, and recalling that m{^) = |||?7^||^2, we 
get 

(55) 5t||w||i2 = -29t/ii9^,m(/ii) - 29i/i2 9^2m(/i2) - 29* Re (77^^,77^2). 

First, using the exponential localization of solitons in space and the fast relative 
phase of the solitons, we estimate an upper bound for 

\dtRe{ri^^,r]^^)\ = \dtUj{r]^-,,e2N+iVa2)\- 

From (l27 l) -( !30l) and ( l49l) . it seems a priori that |9i Re (770-^,770-2)1 is of order ||f |p. 
However, we can pull a factor of ||fi — f2||^^ from the fast oscillating phase 
Q^{vi-v2)-x i^y integrating by parts twice, as in (l23l) - (l26l) in Sect. HI Therefore, 

(56) |aiRe(r/o„r/.2)l<Ce-«""^~'''"- 
Furthermore, (l30l) implies that 

(57) Idtfiid^Mf^i) + dtl^2d^Mf^2)\ < C(/i2 + e"«l''^^-'^^ll + ||i/;||i.)- 
Now, (I55D - (I57j) yield 

(58) IdtMU < C{h^ + e"«ll"^""^ll + llu^lliO, 

for some positive constant C independent of ||t;o|| and h. 
It follows from (!58|l and the Duhamel formula that 

(59) 11^11^2 < C(e^*(/i2 + 11^11^2) + / t/se'=(*-^)e-«ll"^-"^ll). 

Jo 

For times t < C||^;o||^ e ^ (0, 1), we know from (!28|) that 

(60) ||t^i(t)-t;2(t)|| >co|K;o||, 



23 

for some constant Cq > 0. Making the change of variables 

s -^ a(s), 

where 

a{s) := ||ai(s) -a2(s)||, 

and using that 






m^. m\l and dHUD, we have 



t ct 



h \\m\ 

Together with fISU]) . we get the estimate 

(61) \\w\\l2 < C{h^e'' + ^e-^*), 

Iroll 

for some positive constants C and c that are independent of ||fo|| and h. Let 
r := - min(log ||fo||, 2| log/i|) for some a G (0,1). For t < t, (!6T!) imphes 

sup ||w;||i2<C(||i;o|r^+° + /i'(^^")). 
te[o,r) 

n 

7. Proof of Theorem [T] 

We now show that, for ||t;o|| ^ 1 large enough and /i ^ 1 small enough, the 
hypotheses of Propositions [1], [2] and [3] can be simultaneously satisfied. 
Let 

T := sup{t > 0, ip(t) e Ud,n,5 ^ith 6 as in Proposition [T][. 

By continuity of ||w(t)||i2, T > 0. If T < Camin(log ||t;o||7 2| log/i|), then by 
Proposition [3], 

(62) sup \\w{t)\\L2 < C"(||i;o||"^ + h'-"). 
te[o,r) 

Here, C, C appear in Proposition [31 We need 

5 <^ inf m'{fi), 

where 5 appears in Proposition [TJ Consider vq and h satisfying 

C"(||i;o||"^ + /i'-°) < ^<(inf m'(/i))T^ < inf m'(/i). 

Then ||w(T)||i2 < -, andT is not the maximal time unless T = Ca min(log \\vq\\, 2| logh\). 
Then (!62|) yields ([5]). Furthermore, the hypotheses of Proposition [2] are satisfied. 
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Using ([5]) in fl27|) - fl30|) gives the estimates on the evolution of the parameters in 
Theorem [H D 

8. Comments on separating solitons 

We now discuss Remark [2] in Subsect. 11.21 whose hypotheses we assume. 

(1) Suppose that the sohton centers of mass are initially separated by a dis- 
tance d ^ max(— ^, ^L, I loginf^g/g m'(/i)|), such that ||w||i2 = 0{e~^'^) 
for some x > 0, and that ||fi||, ||f2|| = 0(1). Then the analysis above 
(Propositions [H [2] and [3]) holds, except that ( l58l) implies 

Iklli^ < Ce-^'^+^"'^^ < Cd-\ 

for t < d"^, e E (0, 1), from which follows the claim of (1) of Remark [2j 
(2) In the case of escaping solitons, fl58l) in the proof of Proposition [3] implies 

that 

Iklli^ < Ce-^'IWolH-C7"||.„r < c/\\v,\\ 

for t < ||t'o||^) e G (0, 1). Hence the claim (2) of Remark [2] also holds. 

9. Appendix 
Proof of Lemma [U Sect. [3 Explicitly, 

where eaTj^, a = 1, ■ ■ ■ , 2N + 2, are basis vectors of T^^jvi^. For a = 1, ■ ■ ■ , A^ 
and ;3 = iV + 1, ■ ■ ■ , 2iV. 

{e^^.^iep^.) = (-5..Je^^-(^-'^)+*^r/,,(x - a)), -xpe^^-^^-'''>+'^ii,{x - a)) 

V 

= Y(^^M(a;-a),X/3?7^(x-a)) + (9^^?7^(x - a),x^?7^(x - a)). 

It follows from translational invariance of the integral, and positivity and spherical 
symmetry of ^^^(x) that 

{i'nf.ix - a),Xf3r]f,{x - a)) = {ir]i,{x),xpr]^{x)) + (i77^(x), apri^,{x)) = 0, 

and, by integration by parts, 

(5xc'7m(2^ - a),xpr]^{x - a)) = -5„;3(?7^(x),?7^(x)) - {xf3r]^{x - a),d^^r]^{x - a)) 

= -25apm{ij) - (9^.„?7^(x - a),Xf3r]^{x - a)), 

where m{fi) = |||?7^||l2 and 6af3 stands for the Kroenecker delta. Therefore, 

{ea7]^,iep7]^) = -{epr]„,iear]a) = -Sai3m{fi), a = 1,- ■ ■ ,N, P = N + 1,- ■ ■ ,2N. 



25 



Furthermore, 

{d^r]„,id^r]^) = -{d^r]^{x - a), 77^ (x - a)) 

and hence 

where m'{fi) = d^m{n). For a, /3 = 1, ■ ■ ■ , TV, 

where we used translational invariance in the second hne. It follows from spherical 
symmetry of ri^{x) that 

{dx^Vf^, idxfjVf^) = 0, a 7^ /?. 
Furthermore, since 77^ is real, 

(^M'^^m) = {dx^V,^,idx^Vt^) = 0. 
Therefore, 

{eaVa, iefiVa) = 0, a, P = !,■■■ ,N. 
For a = I,--- ,N, 



-•) 5 -^ ' 5 



{ea'n„,ie2N+iria) = ((^^a + dr,Jr]f,{x - a),ri^{x - a)) = 0. 



and 



i 1 

{eaVaJe2N+2Va) = -((^^« + dxJVt^i^ - a),id^r]f,{x - a)) = --t;„m'(/i). 

For a = iV + l,--- ,2iV, 

{ear]a,ie2N+iVa) = -{iXaVfii^ - a),r]^{x - a)) =0, 

and 

{ear]a,ie2N+2VcT) = {iXaVfiix - a),id^r]^{x - a)) = aam'{fi). 
Explicitly, we have 

/ -m{fi)lNxN -^vm'{n)\ 

m(/i) l^vxAf am!{jj) 

m'i^i) 

\ ^v^m'{fi) —a^m'{fi) —m'{fi) / 



^a\%,^Ms 
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where Inxn is the N x N identity matrix, and {■)'^ stands for the transpose of a 
vector in M.^ . One may easily verify that the {2N + 2) x {2N + 2) skew-symmetric 
matrix Q^r given in (fTTj) is invertible if d^j_m{^) > O.D 
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